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ABSTRACT 


Gamma Ray Dose calculations are done using Monte Carlo 
Methods for a spherical shield geometry with a straight 
cylinderical duct. All values obtained are for a disc source 
placed at the center of the shield. All the three major 
gamma ray interaction processess namely Photo electric effect^ 
Pair— Production and Compton scattering are included in the 
simulation of Photon transport. Results are obtained for 
dose points on the duct axis and also for points off the 
duct axis. The latter shows that the scheme of particles 
scoring used in the present work is not adequate and should 


be modified. 
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CHAPTER-1 


INTRODUCTION 

Photon Transport studies are of great importance to 
Nuclear Scientists in shielding analysis and in other allied 
artas such as medicine. In practice^ a few among these shielding 
problems can be handled by mathematical methods. The, then 
obtained analytical solutions or expressions are used to 
estimate the desired quantitites of interest such as dose, 
in a simple way. 

Often Nuclear Engineers fiace problems which cannot be 
dealt by analytical methods because of irreqular shapes in soitrce 
and shield geometries. Shiel d geometry depends on the source 
geometry, the specific, use of radio isotope etc. 

Radiation streaming through ducts is also an important 
aspect in the design of shield. This problem was analysed by 
Tstiro et al (11) analytically. The shield geometry considered 
for their study is an infinite slab shield with finite thickness, 
having a cylinderical duct without any bend in it and the 
source is planar. Apart from research purposes, geonetries with 
openings find extensive use in radiation therapy/non-destructive 
testing. One such geometry has been cho«s^. for the present 
study and the geometry is finite in shape. Description and 
of interest- will be described in chap, 2. 


the problem 
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Though analyt ical solutions ane feasible to some 
extent for duct streaming problems ^ which are, in turn, 
derived with severe assumptions, one cannot have confidence 
over such solutions unless they are verified by seme other 
methods and ultimately by experimental results. The same 
solutions beceme inapplicable when source and shield geometries 
are modified. For solving such radiation transport problems, 
probably Monte Carlo methods are suitable. 

As is well, known that all physical processes including 
the emission of radiation from sources and their subsequent 
transport through material are probabilistic^ in natiare, i,e,, 
one cannot predict with certainity what will happen for 
every individual particle in the process, Monbe Carlo methods 
are used to simulate radiation transport through medi'um 
mathematically, through the use of various probability density 
functions. Density functions are the basic tools of Monte 
Carlo and are discussed in chap, 2, 

For a photon radiation transport problem, two broad 
forms of Monte Carlo method are possible (12), 

(i) Analogue Monte Carlo ; In this method, the 


schematizsibion followed closely resembles with the actual 


3 


physical prcx:ess. 

In this form/ individual particle/histories are 
followed based on mathematical probabilities. Particles are 
created by sampling the governing density function of source 
which describes the behaviour of real particles. Each particle 
is followed until it survivea in the medium of interest/ shield 
for e.g./Ond then a 'new* particle is started from the source. 
Wh^_n large number of particles have been accounted/ the 
acciamulated data properly handled/ will give the quantity of 
interest, so there is a strong analogy _,in this approach^ 
between physical and mathematical particles. 

(ii) Non-Analocfue Monte Carlo; Non— Analogue method 
is more sophisticated to solve a problem without any physical 
intution to real particles' or processes. This method attempts 
to solve the transport equation straightaway. 

Although non-aaalogue method is a versatile approach 
to solve transport problems# it is difficult to form exact 
transport equation for any arbitrary geometry. So more emphasis 
has been shown on analogue Monte Carlo method to investigate 
the chosen problem , 
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The main objective of this study is estimating the 
Gamma ray dose rates at seme points, of interest, for the 
chosen source and shield geometries. In order to facilitate 
the computations, following simplifications are made in 
developing a Monte Carlo code for the problem. 

(i) Gamma, rays interact with matter by more than 
ten distinctive processes. However, only three processes 
are significant to Nuclear Engineering problems and other 
processes can be neglected. These are 1 ) Compton scattering , 
2) Photo-electric effect and 3) Pair -product ion. 

(ii) All scattering interactions aie assumed to be 
isotropic wherever it occirrs in the system. 

(iii) For a history, death can oeexor when the particle 
becomes absorbed, leaves the geometric region of interest or 
loses ita significance owing to other factors (8). 

In the Monte Carlo program, a system built-in function 
random number generator RAN(X) haa bedn used throughout this 
work, v/hich generates random numbers uniformly between 


0 and 1 



CH^PrER<-2 


Mui'iTij I'lill riOj^ S 

This Chapter, deals with the various steps involved 
in the simulation of photon transport. Though there may be 
more than one sampling technique applicable to a perticular 
step^ sampling Compton energy^ e.g,^ here, however, 
relevant sampling methods, say referred by earlier workers* 
are applied for each step of the simulation. Next section 
briefs the description of the problem of interest and the 
necesaar-y steps involved in the simtolation are discussed 
in the subsequent sections, 

2.1 DESCRIPTION OF THE PROBLEM 

Present work concentrates on a problem of calculating 
Gamma ray dose rates at sane points of interest for a source 
and shield geometries cis. shown in Fig 2.1, Source geometry 
is a disc having uniform distribution over it. Shield geometry 
is a solid sphere with a cylinderical duct filled with some 
other material than shield material. Here, main object is 

to compute dose rates along the axis of cylinder and about the 
duct exit also. The latter can be interpreted by computing 
doses at different points on either sidew along the perpendicular 




Fig. 2.2(a) Chosen Co-ordinate System- Spherical Polar Coor- 
dinate (b) Coordinate of a Source Point on the Disc 
Source. 
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to the axis of cylinder , Dose points of interest are shown 
as dots in the same Fig. on the cylinder axis and on its 
perpendicular, also. 

Spherical polar co-ordinate system is chosen to 
represent these, geometries in three dimension and the reference 
axes X/ Y and Z are shown in Fig. 2, 2, a. Disc source lies 
in the X-Z plane and cylinder axis coincides with the reference 
Y axis. Shield radius is measured in m.f.p, units which is the 
value at source energy for the shield material, 

2.2 AN OVERVIEW OF ANALOG MONTE CARLO METHCD 

As mentioned earlier ^ Gamma photons from a source 
incident upon a shield and thereafter^ its consequences 
in the shield can^ however^ be realised mathematically by 
using Analogue. Monte Carlo Techniques, Particles coming from 
a source have to be traced individually i.e.^from birth to 
death or termination. Further divide the shield medium into 
desired number of regions or zones according to the choice 
of interest. The procedure to be followed is shown in 
Figs, 2. 3, a and 2,3 ,b. 

This procedTore may be repeated for a large number of 
histories until the desired statistical acurracy is obtained. 



Fig. 2. 3. a FLOW CHART FOR MONTE CARLO 
MAIN PROGRAM 






NTERM? 


YESV+l 


•’ind the new ene, 
ind direction 
parameters 


RETURN 


Pig. 2.3 .i» PLOW CHART POR SUBROUTINE HISPEE 









10 


All steps giv'»n above are discussed in the next few sections 
in relftventf to the chosen problem, A random process or 
event can be described mathematically by making use of the 
idea of probability density function (pdf ) and cumulative 
distribution function (cdf) which are defined as follows; 

(i) Rcobabilitv Density Function ; p(x) describes the 
relative frequency of occurance of its random variable 
X (12), X can lie anywhere between - oo and + oo, A pdf has 
the following propertiea: 

1) 0 < p(x) for all X 

+00 

and 2) p(x) dx = 1 

(ii) Cumulative Distribution Function .; C(x) gives- the 
probability that a random variable x has a value less than 
or equal to a fixed value. It is defined by 

X 

C(X) = / p(x) dx = prob{X£X) ,..,,(2.1) 

-OD 

cdf has. the following propertiea: 


1 ) 


lim 

X -^ + oo 


C(X) = 1 


2 ) 


lim 

X CD 


C(x) = 0 
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To sample a random value x from a given pdf# first 
select a value C(x) which is •nerely a number# essentially 
from numbers distributed uniformly between 0 and 1# then X 
can be obtained directly from equation ( 2JL ) provided the 
corresponding pdf is a simple one# otherwise some other 
methods are to be followed to saiiple x. 

2,3 SOURCE PARAMETERS 

In order to initiate the problem# one has to start with 
source particles i.e,# photons and their parameters. Because# 
it is essential to know their spatial point of origin and 
direction of motion. A disc source is madeup of number of 
point sources. Since the choice of source parameters for a 
photon is not the same for different source geometries and for 
different source distribution like isotropic, cosine distribution 
etc,# only uniformly distributed disc source and isotropic 
distribution for point source have been assuned’ throughoxit 
this work, 

i) Selection of Space Co-ordinates ; Consider the problem of 


choosing spatial co-ordinates for so-urce particles. The 
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co-ordinates (x, 0 / 2 ) of any point on the disc can be found 

interm a of a selected randtm value R between 0 and radius 

of the disc, and i,- an angle randomly selected between 

- A end + ^ , R and are shown in Fig^ 2,2,b. The 

probability density function describing the distribution 

of aources distributed on a disc as a function of radius 
2 

xs 2 R/ r , where Rj is radius of the disc, and hence, 
the cumulative distribution function from equation (1) 
is (3 ) (s ) 

c(r) = rVrt 

using a random number ' r ' between 0 and 1 instead of C ( R ) 
to select a randan value of R directly from the above 
equation, it gives 


2 , 2 
r = R /Rj 

It follows that 



« « • • ^2 «2 ) 


Frcm the fig. 2.2.b, it can be seen that for a given 
value of R, the source positions are equally distributed 
in azimuthal angle y between - A and +7^ ♦ Then selecting 
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for a random number *r*y it gives 

r = / d^.>/2r = ( r+7^)/27< 

or Y = (2r-l ) (2.3 ) 

Now valUGa of R and ^ are known. So point T(x,z) can be 
located on the disc from the eqxaations x = R Cos vend 
2 = R Sin V . 

ii) Selection of Direction Parameters : After aasigning . 
spatial point of a source particle y or more precisely, after 
locating a point so^Jrce on the disc, next consider the 
problem of assigning direction parameters to source particles. 

A disc source is symmetric w.r. to its axes i.e. it emits 
radiations equally on either sides. Further, a point source emits 
radiations into entire 4 steradians of solid angle. Here 
again ^one can consider only a half of the sphere by taking 
the advantage of symmetry about a plane. So now the problem 
has reduced to a half sphere pointing towaruls +Y direction 
and hence, if a point is chosen on it, that point will give 
the necessary direction to the soxirce particle. Of course 
that point may, however, be expressed internis of direction 
cosines using a polar angle g and an ajsisnvthai angle 0 , 


14 


e and 0 are shown in Fig. 2, 2, a. Selection of 0 and Q can 
be done in the following manner. 

1 ) To select a polan angle ' 9 ' randomly anywhere 

between and 0/ the pdf- interms of cosine of the angle 
'k' - between -1 and +1 - is ^ dk (3 ) and from this pdf, a 
random number r selected for the cdf gives directly/ 



or k = 2r— 1 = ods0 ..*••(2.4) 

2) To choose an azimuthal angle '0' between 0 andA/ 
it is assumed that 0 varies evenly in the interval 0 and 7X . 
SO/ for a randan number 'r'/ choose 0 to be equal to'^r. 

Finally/ direction cosines for the s-ource particle 
in spherical polar co-ordinates are obtained from the 
following equations using the values of 9 and 0 , 

U = Sine Cos0 

Uy = Sine Sir0 

u = Cose . 

z 


(2.5) 
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U^/ Uy and Ug are the direction cosines along XvY and Z 
directions respectively. A subroutine is written to compute 
space and direction co-ordinates of source particles using 
equations (2.2)-(2,5), 

2.4 PATH LENGTH SELECTION 

Thus^ a set of source particle's parameters has been 
selected. Next^ the path length from the source to the point 
of interaction has to be determined. The path length alongwith 
parameters of initial direction cosines defines the point at 
which interaction occurs, But^ two mediums Irave to be 
considered for the geometry. Let duct be filled with air and 
the remaining part of the shield is made of lead. Only 
homogeneous medium has been considered wherever lead or air 
is present in the geometry. So before sampling the path length 
considering the duct^ photons emanating from the soxirce will 
fall into one of the following three cases listed below: 

(i) A photon originating from the source may escape 
through the duct i.e,^ photon without even making any 
interaction with shield comesout of the duct freely passing 
through different mean free paths, as shown in Fig, 2,4, a. 
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This contributes a part to the uncollided or unscattered 
dose at the source energy. One m.f.p, in air at 1.25 MeV 
is approximately 13x10 cm whereas the whole spherical shield 
radius is approximately 9,0cm only and therefore, one can 
neglect the attenuation due to air throughout the particles 
line of flights. 

solving equation of a straight line and that of a 
cylinder will give where as shown in lig. 2, 4. a. 

The co-ordinates of that arbitrary point Q(x' ,y' , 2 ' ) can 
be calculated using the following equations namely, 

X- = X- + U^ki 

y * £= y 4- U 
z* = 2 + 

r 

If y (x' + (y' 5^ + (z* ) is greater than the maximum 

radius of the sphercial shield, th®n that photon will escape 
through the duct, otherwise it will not. 

ii) The situation is shown in Fig. 2.4.b. Photon 
does traverse some distance 'IC 2 ‘-found in the same manner as 
'k^ ' - through the duct and the attenuation in air medium. 
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obviously, can be neglected ^ since Is^ < 3c^* Hence particles 
co-ordinates have to be modified in order to keep track of 
the particle. Substituting ^2 in place of k^ in equation 
(2,6) will yield the latest position of the particle with 
initial energy Eq, from where the particle will travel into 
shield by a distance equal to the selected value of the 
path length. However, direction parameters are retained for 
this purpose. In some other situations i.e,, when is 
greater than duct radius, it may happen to travel some distance 
in air from its origin to the surface of inner cavity sphere. 
This is not shown in the figxure. 

(iii) In this case, photon does traverse no distance 
in air. Of course, this will arise only when both source and 
shield are not separated by any distance from each other 
as shown in Fig, 2,4,c. and ofcourse, the chance of occurance 
of such case is very less unless the geometrical considerations 
are modified. So photon may be allowed to travel into shield 
directly just by sampling the path length. The selection 
procedure for the path length 'S' is as follows. 

The probability of first collision between S and 
S + dg along the line of flight is given by (2) 
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p(s) ds = ^ ds 

where jX is the total lin^r. attenuation coefficient of the 
shield materialvihleh gives the probability of collision 
per iinit length. To sample <rdf is given by 

c.(S) p, dS 

Select a value C(S) using a randwi nunber. r between 0 and 1 
above equation becomes 

X = / 5^'^ ^ds = 1 - 

0 

It follows that 

S = - i In(l-r) 

But 1 - r is distributed in the same manner as r, Ther^fore^ 

S = - i In r (2.7) 

The selected path length 'S' together with 
direction cosines gives the point at which the photon interacts 
with the shield. Interaction may be any of the three events 
namely Cmpton scattering^ Pair- Rcoduct ion and Photo electria 
effect. Sampling techniques applied for these three events 
will be disctussed later. 
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2,5 NEW DIRECTION PARAMETERS 

This section discuss the problem of finding new 
direction of a photon after it encountered a scattering or 
collision event. Consider a Compton scattering event occurs 
at aome point in the shield. Cnee an elastic collision occured 
in a medium / its consequense results change in photons 
direction with losa in its energy. The situation is shown in 
Pig. 2,5 .a. A method of sampling photon energy after collision 
will be discussed in later sections^ . In this section^ 
only particles new direction will be determined for a given 
scattering or deflection angle 0^ and its local azimuthial 
angle as shown in Pig. 2,5^aand 2,5;,b* 0^ and 0^ are 
meas-ured w.r, to particles initial direction. The new direction 
parameters v/ill be very useful to maintain the track of 
a particle from its birth w.ir, to a particular reference axes. 

Having known old values of © and 0 and the values of 

0 and 0 , the new direction parameters U', U ' and U * 

are determined with the help of the following fo;ar equations 

where U ' ^ U ' and U' are the modified direetion cosines 
X y 2 

w.r, to different - axes which are parallel to the 




scattered electron 


Fig. 2. 5 . a Conptcm scattering: an elastic collision between 
a photon and electron. The local orthogonal axes (€.n,C) 
are constructed such that C lies along the original direc- 
tion of the photon's travel. 



(Hcke;-) Particle's 'local' angles of scattering. 8 q is the 
deflection angle, # is the azisiuthal angle. Usually the 
angle • can be assumed to be randomly distributed in the 
range 0 to 2 t. 

Figo 2«5„b Initial and Final Direction of Particle 
in Scattering 
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original reference axes Y and Z rt;spectively. 

Cos©^ = Cos0Co(ffle + slnesinG^Coa 

SinQj^ = yi — Cos©^ Cos©^ 

Cos0j^ = CoaCj^Q - 0 )CdsJ3- Sln(p^ - 0)SiniZf 

and Sin0^ = Sin(j^ - j3)cos0- Cos(j2!^ - iZ!)sinj3 

©^ and are the modified polar and azimuthal angles 
respectively. The above, four equations are taken from Wood 
and Ref, (6), According to Ref, (12), these equations are 
obtained by making use of trignometric relations. Once ©^^ 
and 0 are known, U *# U ' and U ' are readily foiind by 

li X j « 

substituting 0^^, 0 ^ in place of © andjZiin equation (&,5), 
Since the purpose of the above equations is to obtain the 
new direction parameters, much interest is not shown in the 
details, of the derivation of them. The derivation may be 
seen in Ref, (12). In addition to the above method, few other 
methods are also proposed in Refia, (3) and (8), 

2.6 ABCXJT THE DUCT 


With the knowledge of photons new direction parameters 
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and the new energy E ■c sample 

a path length 'S"frcm equation (2^) and calcailate the next 
interaction point using equations given below 

X, = X ' + U' £5* 

Yi = y‘ + s* 

and z, = z’ + U' S' .....(2,8) 

Z- 

The next interaction point may face^ apart from the type 
of interaction, any one of the following three situations 
given below; 

(i) It may lie even outside the sphere or shield 
or so to say, photon escapes from the shield with energy E. 

Of course, this can happen at any time independent of the number 
of collisions. If a photon escapes from the shield then 
ultimately its life history will be terminated. 

(ii) It may lie within the shield. Here it is 
essential to decide the type of event to be aceoionted, 
depending upon that only, photon may either be allowed to 
sxirvive or not to survive. 
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(iii) Photon may encounter a situation as shown in 
Fig. 2,4^.. i.e* photon crossing the duct. Here sampled 
path length 'S'* kept unaltered but it traverse an additional 
distance of 'k" ' in air along that direction while crossing 

1 

the duct. Whv=ther photon crosses the duct or not can be found 
easily frcm the values of 2 ^, So if a photon 

crosses the duct^ then its actual point of interaction will 
be found just by substituting (S' + k" ) instead of S' 
in equation (2,8). 

Further when a photon is allowed to cross the duct 
the point of interaction may or may not lie within the shield. 
In the former case^ photon just crosses the duct and in the 
later, case photon escapes through the duct; ultimately its 
termination. 

2.7 SAMPLING TECHNIQUES 

2,7.1 Ccmpt on Scattering i This effect is the scattering of 
a photon by a free electron which is at rest. In this 
interaction^ photon loses some of its energy from its 
original value to some other value E (E^> E) and is 
deflected by an angle from its original direction of 
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motion (see Fig, 2,3.a), The relationship between E and 
©o ®^P^®ssed interms of Compton wavelengths units namely 
^ = 0,511/ E f where E is the photon energy in MeV^ is 

A— = 1 - Cos©^ 

• o 

where and ^ are the wavelengths before and af ter 
scattering respectively and ( ^ ~^q) gives the increase 
in wavelengths which corresponds to the decrease in energy 
(Eq - E), 


In Monte Carlo calculations^ it is recessary to know 
after a Conpton scattering, the value of A and ©^, Sampling 
from Klein-Nishina (K-N) scattering probability density 
function expressed interms of X = ^ / Aq v/ould suffice this 
requirement , 

Though there are many techniques available to sample 
X frcm the K-N pdf. Wood and Blomqxaist (1) suggest, that the 
sampling technique valid for any incident photon energy and, 
of course, the most widely used also, is that of Khan (7). The 
details of K-N can be seen in Kjef,(12) and the derivation of 




X Value 


fig. flow for Kahn's procedure for sanpling 

the K*N distribution to find the wavelength of a gana pho- 
ton after a Compton scattering. 

Xfi is photon wavelength before scattering (assiaaed known) • 
photon wavelength after scattering and X-R-Xu+j/Xu* 
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Khan sampling technique may be seen in Ref«(l), The procedxore 
for sampling the K-N pdf to find the n&i wavelength ' X ' 
is shown in Big. 2.6. 

Using the selected value of x shown above sampling 
procedure^ Cos©^ can be obtained from the equation 

Cos©^ = 1,0 - >■ X + A 
0*00 

Now the values of > and ©^ are known. Next, the azimuthal 
angle 0^ may be selected in the following manner,. As is 
obvious, 0^ can have any valuea between 0 and 360*^, see 
Big, 2,5,b. This angle may be obtained by sampling at random, 
assuming 0^ is distributed evenly in that range, as 
0^ = 2^r, where r is a random number between 0 and 1. 

2,7.2 Pa ir» Er oduct i on : This is an important reaction for 
photon energies greater than 1,022 MeV, in this reaction'’, 
the incident photon disappears in the field of charged particles 
and an electron - positron pair appeai^with total energy equal 
to that of photon. The energy (E - 1.022 )MeV goes into 
the kinetic energy of electron and positron j^ir. They 
are attenuated in the medium. The electron eventually 


origin*! Y'^pHotBi) 

{£>1.022 lleV) 


0.511 Me 



pair production interaction 
occurs here: positron-electron 
pair created 


tlectron eventually 
absorbed in uteri al 


positron combines with\ 
stray electron and is 
annihilated e ^ e 

V 


0.511 MeV 

Z-annThilation quanta created in 
random, but opposite, directions 
Fig. Scheme of the pair production process. 



The 2-annihilation gamaa photons are aasuned to 
SBterge isotropically* That is, with respect to the ortho- 
gonal axes (x*, y', *'), the arimuthal angle ♦p is uniform- 
ly distributed in the range 0 to 2w, and the cosine of Op 
is uniformly distributed in the range -1 to 1* 

Fig. 2,7.b Direction of Annihilation Gammas Afteir 
Pair- Production 
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absorbed while the positron ccmbines with an electron in the 
attenuating material and in this process two Gamma quanta 
of 0,511 MeV (rest mass of electron) are released. 

In Monte Carlo calculations, for accounting a 
Pair-Production event, two assumptions are made. Firstly, 
annihilation Gammas are released at the position at which 
the Pair-Production occurs. Secondly, the two Gammas are 
released at random in mutually opposite directions. They 
are shown in Fig, 2, 7, a. 

Now, the initial direction of annihilation Gammas- 
has to be sampled. The sampling is done in the following manner# 
Further, it is assumed that the angular distribution of 
Gammas to be isotropic. That is, the cosine of 0^, the polar 
angle, is uniformly distributed in the range -1 to +1 and 
the azimuthal angle 0 is uniformly distributed in the range 

ir 

between 0 and 2;^, and 0 ^ are measured w.r. to parent 
photon's direction and they are shown in Fig, 2,7,b, 

According to Ref, (12), to reduce G,P.U, time, only 
one annihilation Gamma from the pair can be followed; as 
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compensation/ its weight must be mxaltiplied by two* Concept 
of weight will be seen in Chap. 3, The error resulting from 
the above approximation is likely to be Very small because 
of the random orientation with which annihilation quanta are 
released. Finally/ the value of © and 0 are selected from 
the equations . 

Cos© = 2r, - 1,0 

and ~ 

where r-j, and ^2 random numbers. Using these values the 
new direction parameters w.r. to the reference axes, can be 
computed using the eqtiation given in section 2.5. 

2.7,3 Photo Electric Effect : It is an absorptive event and 
in this process, the incident photon is absorbed and 
depending upon its energy, an inner orbital electron (i.e. 
fron k shell) is ejected. The vacancy created is filled by 
a transition of one of the outer shell electrons accompanying 
either a fluorescent emission or an Augar electron. This effect 
is very important for low energy photons in materials with 
high atomic numbers (Z). in this Monte Carlo calculations for 
Photo electric, effect it is assumed that all transitions 
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are acconpanied by the emission of fluorescent lines only, 
Thus^ Augar electrons are completely neglected throughout this 
work. 

The energies corresponding to K and L shells are 
called edge energies because the Photo electric cross aectiion 
^pee^^^ elsewhere decreasing continuously shows a sharp 
discontinuity (edge) at each of them (2), The behaviour at the 
K and L edges are shown in Fig. 2, 8, a. 

The energies and compared with the incident 
photon energy serve to determine the probability of K or L 
shell ejections. So when a photon of enerjgy E > is absorbed 
there is a likely chance of K-shell electron ejection with 
kinetic energy (E - E^). Burthen the vacancy created ia 

ix 

filled by an electron from any one of the outer shells (L^ 

M or N) making a transition accompanying a fluorescent 
emission. The possible K and L fluorescent emission or lines. 

^ ^ shown in Fig. 2.8.b and their 

energies are given in Table 2.1. For elements with ZT'Sl, 

L - fluorescent emission is also possible. 

If the incident energy E is greater than then 
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either: K or L shell ejection is possible. Ofcourse outer shell 
ejection is also aJLlowed. Now, one has to decide which ejection 
is to be considered. So some three nmbers 0-^, 0-^ and 0^, 
which will determine the relative probabilities of K, L and 
outer shell ejections for E define the 

chances of L or outer, shell ejections, for Ej^< E<. Ej^, are 
required. 


The values of 0^, 0j^ and 0-^ are taken from Ref.^S) 
and are shown in Table 2.1. These values are computed 
usdng the Photo electric crosa sections at various edges. 
Further, the yield for a shell i is the total probability 
of emission accompanying electron transition from all outer 
shells to a vacancy in shell i. The individual yields for K 

lines and L - lines are shown in the same fable, a random 

til 

number lying between q_i and implies i shell electron 


ejection and similarly a randan number between 
and Y^/ 0^ implies T** fluorescent emission. Hare also, 
the particle's statistical weight is modified and this concept 
will be discussed in Ghap. 3 . N®ct, the photon with new 
fluorescent energy has to be followed. This, as done before 
of Pair-production, those two angles are sampled in 


in case 
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Table 2.1 Final values of jZ5.^-y 

" ■ I — I, i l A ii 


and 


i 

Shell 

Edge 
Energies 
in MeV 


Y. 

1 

Flu. Ene 
in MeV 

Flu. 

Lines 

1 

0 

0,014699 

0.242 

0,0 

0.0 

0 

2 

L 

0,014699 

1.00 

0,3032 

0.012217 


3. 

K 

0.088004 

2.854 

2.074 

0.074969 

K.1 

4. 

K 

0.088004 

3.957 

3.127 

0.072804 


5 

K 

0.088004 

4.570 

3.712 

0.084777 

Si 

6 

K 

0.088004 

4.715 

3.850 

0.087320 

K,2 


Z = 82 
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the same manner and finally, as usual, the new dix'ection 
parameters can be computed using equations given in section 2,5, 

Although, the method to be followed for a Photo 
electric event is described clearly in Ref a, and (5), 

the derivation and the method of obtaining JZ5^s ' and 'f. s ‘ are 
not explained. 



CHAPTER-S 


RESULTS 

In this Chapter^ result obtained by making use 
of the different sampling techniques that are to be followed 
at various stages of photon transport simulation seen in Chap. 2 
have been presented. All necessary data, are of Israel (10) 
and that of mass absorption coefficient of air or tissue 
are from' Ref. (12), In the development of a computer code 
for the chosen problem# three cases have been considered 
to observe the imparovement of results of one case over 
the other and of course# to attain the closeness of exact 
physical situation also, 

Observatiomhave been made firstly# by forcing all 
interactions to be Compton scattering (cxxn)# secondly# by 
considering only two typ)es of interactions namely Compton 
gQ0^tt.ering and Pair— Production event (com+ppjd).# and thirdly 
by including all the three events viz,# Compton scattering# 
PQi;r;— Production and Photo electric effiect (com+ppd+p)ee)* 

Next section gives a brief idea about the concept of survival 


weights 
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3.,1 CONCEPT OF SURVIVAL WEIGHTS 

In Monte Carlo computations, this idea of statistical 
v/eight has been included to account the absorption effect 
after a scattering. A particle starts from the source with a 
weight of unity. Obviously it has not encountered any 
interaction in the medium so far. The particle weight after 
a collision is the weight before collisioDj multiplied by 
the survival ratio which is the ration of 

macroscopic scattering cross section to the macroscopic total 
attenuation cross section. 

Thus after each scattering event the weight is 
adjusted in the same way to account the effect of absorption 
of photons in the scattering. Here, one can choose R ussian 
Roullette method (8) to decide whether the particle survives 
or not . y 

3.2 COMPUTING DOSE 

As mentioned earlier, the five mfps' shield is 
divided into five regions. Dose conputations have been 
made for each mfp or region at points on the cylinder axis 
and along the circumference of the circle with radius equal 
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to that of spherical shield in the Y— z plane to account for 
the dose pointa along the perpendicular- axis of the duct , 

The right half circle is divided into 18 equal parts w.r. to 
the origin in step of 10 degrees each and hence doses computed 
in each interval corresponds to seme amount of height in the 
perpendicular axis. 

Though results are obtained for all the 18 intervals, 
assuming symmetry one can limit to only either upper half or 
lower half of the semicircle. From the observations, it .is 
found that there is not significant difference in values between 
the two halves. 

Khen a particle crosses a mfp,the associated particle 

2 

current is converted into flux units, say T (no ./cm — sec) 
and then dose rates in mR/hr can be caleulated using the 
relation namely, 

(E) 

Dose Rate = 0.05767 X E X X T. 

where 0.05767 is the constant which converts the unit of 
Pa 

E --- T into dose rate mR/hr units, E is the En'ergy of the 

^ Pa 

on JLn M0V sinci mmmmmm jLS tlno rnsss- o^lDsoxijptxon coof fioiont 

^ 2 

of air or tissue at energy E in cm /gm. 
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3,3 MONTE CARLO OUTPUT 

The three cases mentioned above are given 
below along with the results obtained in each case, 

(i) Here all interactions are forced to be 
Comptons scatterings i,e./ absorptive events cannot occur 
directly. Thus the particle weight is adjusted for by a 
factor of survival ratio after each collision. Results 
obtained are shown in Pigs, 3-, 1 and 3,2, Fig. 3,1 shows the 
distribution of doses rates at each mfp for various angles 
i,e,, along the perimeter of the circle, 

(ii) In thiSy all interactions are forced to be 
either Compton scattering or Pair-Production event. Of course 
if the energy of the photon is less than 1,022 MeV then all 
interactions are ultimately forced to be Compton scatterings. 
If E 1,022 MeV^ then lair-Production is considered to occur 
if r ^O/^Pq where r is a random nunber between 

0 and 1 t otherwise it will not occur and here the survival 
ratio is given by 

Survival ratio = 


As describes in (i) above^ results are shown in 




1 2 3 4 5 

Shield Thickness ( mfps') 


Fig. 3.1 Variation of Dose with Mean Free Path (Com) 


Dose Rate, (mR/hr) 



Fig.3.2 Variation of Dose with Angle ( Com ) 
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Figs. 3.3 and 3,4, Canparision study will be done later, 

(iii) In this case ©hoto Electric, effect is also 
included in the simulation along with other two process 
namely Cempton scattering and Pair- Production. Now one has 
to decide the type of interaction to be allowed to occur 
at the interaction point at sene energy, say E. This can be 
done by comparing a random number 'r' with the ratios 
FpE/y^' Pc/^T >^PPD/^T, ^ between 

and (Say ^ ^ Pc Compton 

scattering will occur. Likewise the choice of the event may 
be decided at any energy E. 

Here also the particle weight after an event is 
modified by the factor survival ratio (JJppp + PepP Ft* 
the case of Photo electric effect, if fluorescent emission is 
allov/ed then the particle will survive with that fluorescent 
energy otherwise its life will be teimiinated. Results obtained 
in the case are showrr in Fig^,3,5 and 3.6, 

3,4 COMPARISION OF MONTS CARLO OUTPUT WITH ANALYTICAL METHODS 

A semi-analytical approach for radiation streaming 
through a cylinderical duct has been made satisfactorily by 
Tsioro et.al. (11). First, this method was applied to the case 



Dose Rate, (mR/hr) 



Shield Thickness (mfps*) 

Fig. 3.3 Variation of Dose with Mean Free Path (Com ♦PPD) 


Dose Rate , (mR/hr) 


10 ^ 


1.25 MeV 
Source 
Histories 



90 80 70 eO 50 40 30 20 X) 

Angle, (deg) 

Fig.3.4 Variation of Dose with Angle (Com* PPD) 




Shield Thickness (mfps') 


Fig. 3.5 Variation of Dose with Mean Free 
Path(Com*PPD^PEE) 



Dose Rate , (mR/hr ) 



Fig.3£ Variation of Dose with Angle (Com^PPD ♦PEE) 
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of Gamma ray streaming after single scattering in a duct for 
a disc source at its entrance. Later ^ it was modified to 
include multiple scatterings also. The dose due to single 
scatterings or singly scattered dose at the duct exit in 
given by 


^ L +00 

.2 / / r f 


= y az/ SL 
° ° Po 


SCA 




r 




where L and Rq are the length and radius of the duct 
respectively 

JPq: initial attenuation coeficient of photon 

attenuation coefficient of photon after scattering 

dg: sampled path length 

dj^: distance between scatcering point and detector 
minus the path across the duct 

Dq: distance between the point of origin and the 
scattering point 

D^s distance between the scattering point and the 
detector 

2 

f s an integral parameter and r = Rq + f / 

: the photon energy after scattering 
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energy absorption coefficlert of air at 

ener^ 

The total dose at any point on the duct axis is the 
sum of unscattered dose and scattered dose ^ and the unscattered 
dose at any point 'L' units away from the disc sovtrce is 
given by (9) 

s 

Vs = “• % ^ <3.2) 

X 0.0567 

where is the flxix and Eq is the initial energy of photon 
( at soxarcG ) , 

Estimation of dose rates by analogue method is as 
follows. Attempts were made first/ to calculate the doses 
along the duct axis and later for points off the axis. As 
such/ no point detector calaulation/ say a particiHar point in 
the shield/ can be done straight away Using this method 
but the normal procedure is to consider some distance +a.d 
about that point along the boundary so that the sxm of the 
individual contributions of each particle's dose within that 
total 2Zid length will be considered/ effectively/ as the 
dose about that point, Cfcourse/ the same proceduf-e has to 
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be followed for all points of interest. Finally one can 
get a histogram, say dose vs shield thickness, and from that, 
distribution of doses along the points of interest may easily 
be found. 


Now, the question arise, is, how to choose the 
interval .ad for a given problem”?. As such, there is no specific 
rules or procedvures regarding this. Perhaps some experienced 
workers might guess the procedure to be followed for a problem 
of interest by intution. But in this work, few approaches 
have been made and the results obtained are shown in Figs, 3.1 
to 3,6. 

In the first approach 0,5 and 1.0 were considered for 
/id values i.fc. Ad along Z axis (refer Fig, 2,1) in two 
trials. It was found unsatisfactory for even points on the duct 
axis because the values of total dose cemputed minus the 
unscattered dose computed from equation 3,2 were much higher 
than the unscattered doses (values are not shown here). 

Further some results were obtained for points off 
the axis also using the same approach. Here, the same 2/s.d 
length cannot be considered for every 2/>d, the arc length 
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(boundary) is not coA£±.ant, So it is better to consider equal 
angles/ i.e. divide each mfp into eighteen 10 degree intervals, 
instead of equal heights. Though this helped to remove the 
difficulty faced due to the unequal arc lengths, again, the 
magnitudes of doses were high. 

In the noxt approach, in addition to^iid value, i*d 
i.e, along X axis also included to observe the effect on 
dose values. ^d_ = 0.5 and A d„ = 0.5, showed erratic results. 

Z X 

later when increased from 0,5 to 1.0 (equal to the 

radius, of the duct), satisfactory resulta have been obtained 
as shown in Table 3,1 for dose points on the duct axis. 

Assuming that the above fixed values of end 

dd are valid for points off the axis also, observation 
2 

have been maiie for various angles and the results are shown 
in Figs. 3.1, 3.3 and 3.5, For e.g.. Fig, 3.1 shows the 
distribution of dose rates for various shield thicknesses 
at each angiilar intervals, while Figs. 3,2, 3,4 and 3,6 
show the distribution of dose rates for various points on the 
surface of the shield in each mf ps ' . 

The computed multiple scattered doses on the axis 
are shown in Table 3.1 for the three different cases, mentioned 


Table 3,1 Comparis^on of Monte Carlo Dose Rates with 


o o o o 


u 
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in section 3.3^ considered for obtaining these results. These 
values are compared with the singly scattered doses computed 
using equation (3 ,1), 

For water at 1.00 MeV photons (survival ratio 1.0 
and for 1/Rq =10; the ratio of maximum shield thickness 
to duct radius), it i s seen that (9) the ratio between 
scattered and unscattered dose is about 40%, 3D,for lead, 
at 1.25 MeV (sxirvival ratio approximately 0.8) and for V'Rq 9*0, 
it is expected that approximately the same 40% will be the 
scattered dose. Table 3.2 shows not much difference between 
the Monte Carlo output and the assumed values. 

l\s. is known that monoenergetic photons incident 
upon a shield tmergea from the shield with a continuous energy 
say from 0 to Eq ( the maximum energy Eq). So emphasis has 
been shown to count the number of particles falling in an 
energy interval, say in step of 0.1 MeV from 0.0 to 1,25 Mev, 
while crossing a region or boundary of interest. Fuirther 
this observation v;ill put light on the difference between 
one case and the other (see section 3.3). The observed 
number of particles for different energies at each mfp are 
shown in Pigs. 3.7 and 3,8. 
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3.5 DISCUSSIONS 


In the Figc. 3.1, 3.3 and 3,5 curves are not drawn 
for angles 20°, 30° and 40° because the computed values for 
these; angles lie close to each other. The computed value of 
dose at- the fifth mfp at angle 10° is 0.0 in all the three 
cases considered and this point is not shown in all the three 
figures mentioned above. In practice, there is seme contribution 
from a disc source to that point and hence this gives rise 
to doubts about the choice of the csd^. and values chosen. 

A suggested approach is given in Chap, 4, 


Further, number of particles escaping through the 
duct (i.e. part of the uncollided flux) is found to be 
approximately 620 in all the three cases considered for 10 
source particles. Since the source energy is only 1,25 MeV 
there is not much difference between curves in Figs. 3.1 and 

3.3, and. Figs. 3.2 and 3.4 because Photo electric effect 

^ noi- included. But the same 

which accounts for absorptiOTi xs- ri a-” 

. . o c and 3-6 which include 

figures compared with Figs, 3.5 

a ■i-p-Pt-:iY~^nce in the dose values. 
Photo electric effect show some dif 

the nxmber of particles having 


As mentioned earlier , 



No- of Particles 



Rg. 3.7 Variation of No. of Particles with Mean 
Free Path. 


No. of Particles 



Fig. 3.8 


Variation of No- of Particles with Mean 
Free Path. 
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Table 3.2 Comparis-.on of Monte Carlo Dose Rates 



with 

modified 

Dose Rates. (40% 

above unscattered) 

MPP 

Distance 

In cm 

Dose* 

mR/hr 

Dose M.C. 
case (i) 
mR/hr 

Dose M.C, 
case (ii) 
mR/hr 

Dose M.C. 
case (iii ) 
mR/hr 

1 

2,547 


9,392v- 

9.6450 

9.417 

2 

4.0931 

3,9203 

3.770 

3.9740 

3.701 

3 

5.6396 

2.0933 

2.050 

2.1910 

1,986 

4 

7.1861 

1,2967 

1,351 

1.3470 

1.233 

5 

8.7327 

0.8809 

0.892 

0.9033 

0.875’ 


* Unscattered Dose +40% Unscatter t:d Dose 
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particular energy as a function of distance, in mfps* are 
shown in Pigs, 3,7 and 3,8, Pour sets of exurves fear energies 
0,3, 0*5, 0,6 and 0,9 MeV are shown in those figures. Here 
again, due to absorption in all the four sets^; the inclusion 
of Photo electric effect shows lower values. The amount 
decreased from that of the other two cases is about 30% to 
65%, At 0.5 MeV,the contribution due to Pair— Production is 
more compared to other energies because of the annihilation 
gammas 



CHAPTER- 4 


CONCLUSIONS AND RECOMMENDATIONS 


4.1 CONCLUSIONS 

Tho following conclusions can bt drawn based upon 
the results of the present study. 

1 , Results have been obtained by using Monte Carlo 
methods for the investigation of photon transport 
in a chosen shield geometry shown in Fig, 2,1. 

Computed Gamma ray dose rates along the axis of the 
duct are shown in Table 3,1, Those values have good 
agreement with the results obtained by Tsuro et.’aii(ll), 

2, No conparision has- been made for the doses computed 

at points off the duct axis because neither analytical 
solutions nor experimental results are available, 

3 , In the simulation of photon transport^ all the three 
events Photo electric effect, mir-Production and 
Compton scattering are included and the effects of 
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absorption due to Photo electric effect and E^ir- 
Production are found clearly from Figs. 3,7 and 3,8. 

4, From the values obtained for points, in the source 

plane, (refer section 3,5), one can conclude that 
the values of a dg and A. are not reliable, 

4.2 RECOMMENDATIONS 

1, For the dose calculations, the values of /^d^ andAd^^ 
may be changed or the whole scheme of particles 
scoring should be modified suitably to calculate the 
flex accurately. This can be done, for e.g,, as 
follows. Divide the perimeter of the semi circle 
into 18 intervals, try to make a circle about each 
such interval so that the area under each circle 
remains constant and hence, the number of particles 
crossing each circle can be calculated, 

2, This investigation can be carried out by including 
seme biasing techniques for path length selection 
and depending upon the interest, directional biasing 
may also be included. 
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3 . 


Simulation approach can bo extended to shields made 
of heterogeneous materials or to multilayered shields. 
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